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The setting of Mossbauer nuclei embedded in thin-film cavities has facilitated an aspiring plat¬ 
form for x-ray quantum optics as shown in several recent experiments. Here, we generalize the 
theoretical model of this platform that we developed earlier [Phys. Rev. A 88, 043828 (2013)]. The 
theory description is extended to cover multiple nuclear ensembles and multiple modes in the cavity. 
While the extensions separately do not lead to qualitatively new features, their combination gives 
rise to cooperative effects between the different nuclear ensembles and distinct spectral signatures 
in the observables. A related experiment by Rohlsberger et al. [Nature 482, 199 (2012)] is success¬ 
fully modeled, the scalings derived with semiclassical methods are reproduced, and a microscopic 
understanding of the setting is obtained with our quantum mechanical description. 

PACS numbers: 42.50.Pq, 42.50.Nn, 42.50.Gy, 76.80.+y 


I. INTRODUCTION 

With the advent of novel light sources, the emerging 
field of x-ray quantum optics has gained considerable 
momentum, both on the experimental and theoretical 
side pfl. While ideas based on quantum coherence and 
interference could in principle be realized with inner-shell 
electrons, solid state targets obeying the Mossbauer ef¬ 
fect [2] have sparked interest in many recent works US]. 

A particularly interesting setting in which quantum 
optics with Mossbauer nuclei can be realized is specifi¬ 
cally engineered planar x-ray cavities. Embedding a thin 
layer of resonant nuclei in such cavities has facilitated the 
observation of a number of phenomena, such as the coop¬ 
erative Lamb shift and single-photon superradiance E3, 
Fano line shape control and interferometric phase mea¬ 
surements HU, magnetically controlled reflection spectra 
modified by spontaneously generated coherences |12] and 
group velocity control of x-ray pulses m- 

However, there is an additional cavity configuration 
which has sparked interest recently. In a setting with 
two particularly placed ensembles of Mossbauer nuclei 
in the cavity, the iron isotope 57 Fe with its transition 
at 14.4 keV, it was possible to observe a reflection spec¬ 
trum with a deep interference minimum in the center due 
to the phenomenon of electromagnetically induced trans¬ 
parency (EIT) MM- This is a remarkable result, since 
typically two coherent driving fields are required for this 
effect to emerge. In Ref. ns], however, the EIT exper¬ 
iment was established in a thin-film cavity with only a 
single excitation from a synchrotron beam, whereas the 
second field was intrinsically provided by intracavity cou¬ 
plings between the two 57 Fe layers. For the modeling of 
the experiment, different semiclassical approaches can be 
employed nanBi, while a consistent description based on 
quantum optics as usually desired for EIT is still lacking. 

A first quantum optical model for the light-matter in¬ 
teraction in the cavity was developed in Ref. m- How¬ 
ever, it does not yet cover cavity settings with multiple 
resonant layers, and hence it is also not yet capable of 
describing the EIT experiment. But motivated by the ex¬ 


pected significance of multilayer configurations, it would 
be highly desirable to also have a microscopic theory at 
hand, which allows for a deeper understanding. Top¬ 
ics of interest include the nature of the above-mentioned 
intrinsic cavity-mediated coupling between the different 
resonant layers, and perspectives on how multilayer cav¬ 
ities can be specifically engineered. We note that other 
approaches based on scattering theory were employed to 
model the two-layer layout |2(T| : however to our knowl¬ 
edge they remained unsuccessful in providing a quanti¬ 
tative description of the EIT experiment. 

In this work, we generalize the quantum optical the¬ 
ory from Ref. m with the aim to describe the single¬ 
photon EIT experiment m and related settings. To this 
end, we extend the description to include multiple cav¬ 
ity modes as well as multiple layers. The extension to 
multiple modes allows us to accurately describe the cav¬ 
ity reflection in the absence of resonant nuclei. How¬ 
ever, the general shape of the nuclear contribution to the 
measured signal turns out to be unchanged when con¬ 
sidering the two extensions separately. In this case, the 
line shape predicted in the absence of a magnetic hyper- 
fine splitting is a Lorentz profile, in which only the co¬ 
efficients are modified due to the additional elements in 
the theory. However, if multiple cavity modes and mul¬ 
tiple layers are considered simultaneously, a new class 
of nuclear reflection spectra is obtained. In particular, 
restricting the analysis to two resonant layers, the EIT- 
like spectrum observed in Ref. m is reobtained and the 
predicted scalings are in accordance with previous semi¬ 
classical calculations. Due to the microscopic ansatz of 
our theoretical model, we can provide a full quantum 
interpretation of the system. A good agreement to nu¬ 
merical results obtained from semiclassical descriptions is 
found over a broad parameter range. Furthermore, our 
extended model opens up avenues to engineer a broader 
set of effective level schemes at x-ray energies by com¬ 
bining the advanced possibilities of multiple modes and 
layers together with magnetic hyperfine splitting. 

This paper is structured as follows: In Sec. |TT] we re¬ 
capitulate the cavity system and the basic model which 
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FIG. 1. (Color online) Schematic of the considered setup. A 
thin-film cavity is probed by hard x rays in grazing incidence. 
The light-matter interaction with the resonant nuclei in the 
center of the cavity, possibly under the influence of a magnetic 
field, modifies the observed reflected signal. 
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was developed for its description in Ref. m■ After this, 
we generalize the model to include multiple layers and 
multiple modes. In Secs. ||V] and [VJ we analyze the ef¬ 
fects of the two extensions separately and in Sec. |VI| the 
consequence of both extensions applied simultaneously 
are discussed. Finally, the general model is applied to 
describe the EIT setting from Ref. |16] . 


FIG. 2. (Color online) Field distribution in the two cavities 
defined in Tab.|I]is shown. The incidence angle (9 ~ 3.5 mrad) 
is chosen such that the third guided mode is driven. The cav¬ 
ity field is normalized to an input field with unity intensity. 
Solid red lines take into account only electronic scattering in 
the cavity, which is realized at x-ray frequencies off-resonant 
to the nuclear transition in 57 Fe, dashed blue curves assume 
resonant driving of the transition. The shaded areas (red, 
gray, blue) indicate the positions of the different cavity ma¬ 
terials (Pt, C, Fe). 


II. RECAP OF THE BASIC MODEL 

The setup which is considered in this work is a thin- 
film cavity with embedded Mossbauer nuclei, probed by 
hard x rays in grazing incidence as visualized in Fig. |T] 
Such a cavity is formed by a stack of different materials. 
At its boundaries materials with a high electron density, 
such as platinum or palladium, act as mirrors, while the 
material in the center, e.g. carbon, has a low electron 
density and provides a guiding layer for the x rays. For 
certain incident angles in the mrad range, the x rays can 
resonantly excite a guided cavity mode and propagate 
inside the cavity, rendering the structure a waveguide-like 
system. By embedding Mossbauer nuclei in the center 
of the waveguide, a near-resonant interaction of the x- 
ray light with the transitions of the resonant nuclei is 
achieved. The reflected signal forms the main observable, 
and its spectral shape is crucially influenced by the light- 
matter interaction in the cavity. As proven in a number 
of recent experiments fTUHT3l fl6] , this setup constitutes 
an auspicious platform for the exploration of quantum 
optical phenomena in the x-ray regime. 

A quantum optical model for the description of the 
x-ray light-matter interaction in these thin-film cavities 
was introduced in Ref. HE however, it is limited to a 
small subset of possible cavity layouts. In the following, 
we will briefly present the existing theoretical approach, 
before we then continue to generalize the model in order 
to cover more elaborate scenarios. 


A. Cavity 

As already mentioned above, guided cavity modes can 
be excited for certain resonance angles 9q. These angles 
depend on the cavity layout, such as the materials and 
the layer thicknesses, and can be determined by comput¬ 
ing the angular-dependent reflectance in the absence of 
resonant nuclei. At certain positions of this curve, the 
reflection is strongly suppressed, indicating the presence 
of a guided cavity mode. The reason for this suppression 
is the destructive interference between the reflection di¬ 
rectly at the cavity surface and the reflection of the light 
which entered the cavity mode. 

Exemplary field distributions for two particular cavi¬ 
ties, which will be analyzed later, are shown in Fig. [2] 
The incidence angle is chosen such that the third guided 
mode is driven, which is reflected in the three antin¬ 
odes of the field intensity inside the cavity. The external 
cavity field stems from the interference of the incident 
(“exp {ik z zy) and reflected beam (“Rexp (— ik z z)”). 
Hence, the suppression in the reflectance R <C 1, which is 
characteristic for the guided modes, is visible as a small 
modulation of the external field. In contrast, if the x- 
ray frequency matches the transition of the nuclei in the 
cavity (dashed lines in Fig. [2]), the relative strength of 
the reflected light becomes larger as it now consists not 
only of electronic, but also nuclear scattering contribu¬ 
tions which adds up to the observed signal. Note that 
the external field in Fig. [2] is shown close to the cavity 
surface and therefore it an interference pattern of the in¬ 
cident and the reflected light appears. In the far field, 
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TABLE I. Parameters for the two-layer cavities analyzed in 
this work. The geometry defined in the first column corre¬ 
sponds to a node-antinode configuration in which an EIT- 
like spectrum is expected, the second parameter set defines a 
antinode-node configuration in which no EIT-like spectrum 
was observed M- 


Material 

Thickness [nm] 

Thickness [nm] 

Pt 

3 

3 

C 

10.5 

20 

57 Fe 

3 

3 

C 

6.5 

6.5 

57 Fe 

3 

3 

C 

21 

11.5 

Pt 

10 

10 


however, these two contributions are easily distinguished 
due to the different propagation direction, as visualized 
in the schematic in Fig. |T] 

Since the different cavity modes are well separated in 
their resonant incident angles by several mrad, only at 
most one mode is usually driven near its resonance. Nev¬ 
ertheless, to take into account the two polarization states 
of the x rays perpendicular to the propagation direction, 
two cavity modes a\ and ci 2 are included in the theoretical 
description. In the rotating frame of the external driving 
field with frequency w, the Hamiltonian characterizing 
the free evolution of the cavity modes with photon anni¬ 
hilation and creation operators a and a) as well as their 
coupling to the external field a- m is given by (h = 1 used 
here and in the following) [T9] 


Hm = Ac-ai^ai + Ac-a^ &2 

+ W (al-a in )a in a[ - (a* n -ai) a* n ai 
+ W2k r (a 2 -d in )ai„a£ - (a* n -a 2 )af n a 2 . 


( 1 ) 


Here, the expressions in the round brackets denote scalar 
products of the polarization directions between the inci¬ 
dent radiation (di n ) and the cavity modes ( 61 , 02 ). The 
mismatch between the frequencies of the cavity mode and 
the external field is denoted by the cavity detuning 


A c(9) = wc — u> = w 


\]cos ( 9 ) 2 + sin (9q) 2 — 1 


( 2 ) 


for arbitrary operators 0 + and O . Then, the photon 
loss can be described via 

£m[p\ = - KC[p,a\,ai\ - nC[p,a\,a 2 \ ■ (4) 


B. Nuclei 

Next, we include the resonant nuclei to the descrip¬ 
tion. In this work, to be specific, we refer to the com¬ 
monly used Mossbauer isotope 57 Fe with its transition at 
Wo = 14.4 keV and a line width of 7 = 4.7 neV. To en¬ 
compass the general case, we allow for a magnetic hyper- 
fine splitting of the nuclear resonance. Then, the ground 
state splits up into two states separated by the energy 5 g 
and four excited states with energy spacing S e between 
adjacent states. This leads to six Ml allowed transitions 
in the nucleus [2T j. which are summarized in Tab. |ll| We 
note that the energy splitting of the ground states is in 
the neV range, such that both states are evenly popu¬ 
lated at room temperature according to the Boltzmann 
factor exp (Sg/ksT). 

In a suitable interaction picture, the nuclear dynamics 
is characterized by the Hamiltonian 


H n = H 0 + H c . 


( 5 ) 


The diagonal part 


N r 2 

H o = Yl I) \9j n) )(9j n) \ 

n =1 '-j =1 

+ E («-§>-A) 

j= 1 


( 6 ) 


contains the energy detuning A = w — ujq and the energy 
splitting of the states due to the magnetic hyperfine in¬ 
teraction for all N atoms. It is important to note that 
in contrast to the cavity detuning Ac, the detuning A 
does not depend on the incidence angle 9, but only on 
the frequency of the externally applied x-ray field. In 
Eq. ([ 6 ]), the index n sums over all N nuclei, and j sums 
over the two ground and four excited hyperfine states, 
respectively. The coupling of the six transitions p (see 
Tab. 0 to the two polarization modes are given as 


We emphasize that the cavity detuning can be controlled 
with the incidence angle 9 and does not depend on the 
frequency of the driving field for all practical purposes. 
The reason for this is the special grazing incidence geom¬ 
etry, in which the cavity is probed, and that the nuclear 
resonance typically is orders of magnitude more narrow 
than the cavity line width. 

Furthermore, incoherent effects such as the photon loss 
of the cavity modes have to be included in the model. 
This can be done via a description in terms of the density 
matrix p and Lindblad operators C[p\. We define 

C\p, 0+,0~] = ( 0 + 0~p + pO+O- - 2 C>-pO + ) (3) 


N 6 


n —1 /i =1 


+ (<-fi 2 )ff£ n) sS «2 + H.C.1 . (7) 


~ • • (n) 

Here, d^ denotes the normalized dipole moment, S y ^_ 

[S 1 ^] is the nuclear lowering [raising] operator of the 
transition p for atom n, and the coupling coefficient 
= gcg,e l ^ is composed of a constant g , the Clebsch- 
Gordan coefficient of the respective transition and a 
phase taking into account the position of the atom n. 
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TABLE II. Overview of the Ml allowed transitions in the 
57 Fe nucleus with transition index p. Shown are the involved 
states, the transition energy A E relative to the energy at 
vanishing magnetization u>o, the Clebsch-Gordan coefficient 
(CG) c M and the polarization type. Linear polarization is 
denoted by n°, right (left) circular polarization as <r + (cr _ ). 


9 

Transition 

A E 

C-G 

Polarization 

1 

Iffi) +* |ei) 

-5J2 - 3/2c5 e 

1 

<J~ 

2 

Iffl) ^ l e 2> 

-5 g /2 - l/28 e 

^2/3 

7T° 

3 

Iffi) ^ |e 3 > 

-S g /2 + 1/2 5 e 

yi/3 

(7 + 

4 

|S 2 > |e 2 ) 

8 g /2 - l/2S e 

yi/3 

a~ 

5 

1 92 ) ** |e 3 ) 

8 g /2 + l/2S e 

vw 

7T° 

6 

Iff 2 > ^ |e 4 > 

Sg/2 + 3/2<5 e 

1 

<7 + 


Additionally, the spontaneous decay of the excited nu¬ 
clear states is taken into account via 

N 

£se[p] = E4e } [p], (8) 

n—1 

4bM = ■ P) 

A*=1 


C. Input-output relation and observables 


In order to relate the internal operators in the cav¬ 
ity to externally accessible quantities, the input-output 
relations are employed (23]. The output field a ou t, also 
visualized in Fig. |TJ is given by 

t^out — CL'm (a ou t-«in) 

+ \/2K^[(a* ut -ai)ai + (a* ut -a 2 )a 2 ] . (10) 


With this operator at hand, the reflection coefficient 
reads 


( a out) 
ttin 


(ii) 


Note that in a typical experiment, the reflectance |R| 2 is 
measured. 


D. Full model 

The expressions given above form the building blocks 
of the general model developed in Ref. [15]. The full 
master equation reads 

~^p = — i[Hm + H n , p] + Cm[p\ + £se[p] • (12) 

In principle, the dynamics of the system could be solved 
this way. However, due to the huge Hilbert space con¬ 
nected with the N atoms and two cavity modes, this 


task is challenging. Hence, in Ref. m two approxi¬ 
mations well justified at present experimental conditions 
were performed. First, the cavity modes ai and ct 2 were 
adiabatically eliminated. This is possible since the cavity 
modes have a low quality factor Q im which is known as 
the bad-cavity regime [23]. Hence, the time scale 1/re, on 
which the mode dynamics equilibrates, is very short com¬ 
pared to the nuclear time scale and therefore the cavity 
modes can be considered as stationary. Second, it is pos¬ 
sible to restrict the analysis to the subspace of up to one 
excitation in the system, since experiments performed at 
current synchrotron radiation sources provide on average 
less than one resonant photon per pulse mmm - This 
simplifies the master equation considerably and compact 
analytic expressions for the observables can be obtained. 
Below we will exploit the same approaches to also sim¬ 
plify the extended model which will then include different 
resonant layers and multiple modes. 


III. GENERALIZATION TO MULTIPLE 
LAYERS AND MULTIPLE MODES 

We will now extend the basic model introduced in the 
last section by explicitly including multiple layers of res¬ 
onant nuclei as well as more than a single cavity mode 
in theory. As already discussed above, in typical experi¬ 
ments at most one cavity mode can be driven resonantly 
at a time. This is due to their large angular separation, 
which significantly exceeds the beam divergence at mod¬ 
ern x-ray sources. The angular separation in turn leads to 
large cavity detunings Ac for nonresonant modes. Nev¬ 
ertheless, as we will show below, the additional modes 
can sometimes be of importance, since the nuclei can in 
principle scatter into them, or if the reflectance is con¬ 
sidered over a broad range of incidence angles. Also the 
inclusion of multiple layers is a highly desirable goal, as 
motivated by observation of EIT in Ref. m- 

Clearly, the coefficients in the master equation, such 
as decay rates re or coupling constants g , differ for each 
mode and each layer and ought to be marked with an in¬ 
dex for the respective element. In an attempt to reduce 
confusion, we stick to the following notation: An atomic 
index is denoted by an upper index n in brackets. Lower 
indices p indicate a transition, as listed in Tab. [TT| The 
different cavity modes are distinguished by an upper in¬ 
dex \j\ in squared brackets. A curly bracket {?} indicates 
that the respective quantity is related to layer l. This 
notation is summarized in Tab. m 

We start by revisiting the internal and external elec¬ 
tromagnetic field. We consider a single incident field a- lni 
which impinges onto the cavity surface under the graz¬ 
ing angle 9 , and an outgoing field a out , emitted at the 
respective reflection angle 7r — 9. Compared to the ini¬ 
tial analysis in Sec. [IT] the input held does not only drive 
one cavity mode a, but multiple modes a^. At the same 
time, the output held is driven by these modes and, nat¬ 
urally, also the resonant nuclei will interact with the dif- 









5 


TABLE III. Index notation used throughout this work. 


Index 

Quantity 

n, m 

nucleus n, m € {1,..., N} 

fl,U 

nuclear transition g, v £ {1,. . . , 6} 


nuclear layer 

L?1 

cavity mode 


ferent cavity field modes. This means, we also have to 
distinguish the coupling coefficients and decay rates for 
each cavity mode. 


Generalizing Eq. (10) from the original theory, we write 


for the input-output relation 


_ / ~ * 

^out — ^in (a oU 




2k 


[j] 


Li] (ft* 

n 1 u o 


Aj] 




(“out '“2 1 ) 


(13) 


and the Hamiltonian describing the dynamics of the 
modes given in Eq. |l]) becomes 


nuclei is present in the cavity. Here, the nuclei close to 
the two layer boundaries might be exposed to strongly 
differing field strengths and the respective coupling coef¬ 
ficients become spatially dependent. 

Both cases can be modeled by introducing several en¬ 
sembles of nuclei. The atoms in each ensemble are sit¬ 
uated at the same depth of the cavity and hence couple 
to the modes with a common coefficient. We denote this 
coupling parameter between the nuclei in the layer l and 
the cavity modes j by g LfiW. The coupling coefficient of 
the transition g in single atom n located in the layer l n 
then reads 


9, 


(n)b'] _ l c 


= 9 




(15) 


The number of nuclei in each layer is and the to¬ 

tal number of resonant nuclei is N = ^./V^. In this 
formulation, the coupling Hamiltonian from Eq. ([7j) is 
generalized to 


(<W) 9? )M sS“S i] 


n,Li,3 


H.C. 


(16) 


Hm = J2 Al3] ^ 


a in a 


CL i CL- — l - CL<2 CL<2 


i Y\ 2K 


AJ] 

'R 


1 ] (“i Jl -“in) (a* n -a^ ] ) 

(“ 2 1 -“in) -a*na [ 2 


+ a in af(^* 


(14) 


In a similar fashion, the couplings with nuclei are mod¬ 
ified to include the sum over all modes j and the in¬ 
teraction Hamiltonian given in Eq. ([ 7 ]) and describing a 
transition g of an atom n is extended accordingly. While 
the coupling to the single cavity mode was denoted by 
gl'"' 1 before, the coefficients related to a general cavity 
mode j are now named g M 1 . They can be decomposed 
into = g^Cf l e t< i >M , where is the Clebsch-Gordan 

coefficient of the transition g, </>(") accounts for a poten¬ 
tial phase imprinted on the nucleus by the field due to 
the atomic position, and g denotes a universal coupling 
constant between mode j and all nuclei and transitions. 
Note that this factorization is possible in this way only 
because we assumed a single thin layer of resonant nuclei 
in this analysis. 

However, as soon as we consider multiple layers, the as¬ 
sumption of uniform coupling strengths g ^ for all atoms 
is clearly no longer justified. For example, for a cav¬ 
ity without resonant nuclei, roughly an intensity profile 
with sin 2 shape can be expected along the cavity for the 
guided modes. Different layers at different positions will 
thus experience different field strengths and the coupling 
coefficient g to the cavity mode cannot be considered as 
a constant anymore. Also, we want to emphasize that 
the same argument holds if a very thick layer of resonant 


The diagonal part Hq containing the energy shifts of the 
states and the detuning A is unaffected by our extension 
of the model. 

Next to the Hamiltonian dynamics, also the incoherent 
part capturing the mode decays need to be extended ac¬ 
cordingly. The Lindblad operator describing the photon 
loss in the cavity modes, see Eq. 0 , becomes 

£m[p] = ~Y nb] c iP> a i ] ^ a i ] ] 

3 

- Y £[p, a b] \ a 2 ] l , (17) 

J 

whereas the spontaneous emission contribution of the nu¬ 
clei remains the same. 


A. Effective Master equation 


In a next step we simplify the master equation by ap¬ 
plying the same approximations as in the case of the orig¬ 


inal model, which were described already in Sec. IID 


First, we perform the adiabatic elimination of the cav¬ 
ity modes. In contrast to the basic model, we do not 
eliminate the two modes ai and 02 for the two polar¬ 
ization directions only, but a total of 2 j modes. How¬ 
ever, since the different modes are not directly mutually 
coupled, they can be eliminated independently and their 
contributions to the effective master equation sum up. 

From the Heisenberg equation of motion for the cavity 
mode operators 


dt 


= i[H M + H 0 + H c , «[ 2] ] - K [i] a [j] 


(18) 
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we find the stationary solutions 


aP = 


1 


«b'] + iA^ 


2K$a in (a/ i -din) 


IjV 




n,fi 


(19) 


where i = 1,2 indicates the two perpendicular cavity 
mode polarizations. 

Inserting these operators in the full model, we obtain 
the effective master equation for the nuclei 


^ P — * [Hefi i p\ T ^eff \p\ 


( 20 ) 


with the effective Hamiltonian and the Lindblad terms 


H e s = H 0 + Hq + H ls , 

£eff[/?] = £se[p] + £cav[p] • 


( 21 ) 

( 22 ) 


In the same notation as in Ref. m, the individual com¬ 
ponents of these equations are found as 

= £ E (nWjWW) s<" + > + h.c. , 

3 

(23) 

»LS = EE «-‘x-4) E (W“9f”>“ r ) 


n,m 

n(n) Q (m) 


(24) 


cw=EE (<-‘4 a) E (-^^“sf”®') 

3 

(25) 


with the coefficients 


= 


2 K 2 ^in 


«b‘l + iA^Q 


^=Im 


CF = Re 


Atbl + iA 1 /} 


.bl 
V; , 


(26) 

(27) 

(28) 


and the outer product lj_ = didj + Note that 

this completeness relation only refers to the two possible 
mode polarizations, such that no sum over the different 
modes \j\ is required. Moreover, the adiabatic elimina¬ 
tion also affects the input-output relation from Eq. (10) 
and thus the observable in Eq. We obtain for the 

reflection coefficient 


(29) 


with the cavity contribution Rc and the nuclear part of 
the reflectance Rjy. The two reflection coefficients are 
given by 


R c = “1 + E 


2k 


bl 


«bl + iA [ / ] 


(“out'“in) i 


(30) 


^ = E 


X (S™) - 


2 nWgPW 


cbl 


iA [ /J 


(“out' 1-t, 


■d, 


(31) 


In a second approximation, we restrict the dynamics 
of the system to the subspace of up to one excitation in 
the system. As mentioned above, the reduction to the 
linear regime is well justified for experiments performed 
at current synchrotron radiation sources. In the initial 
stage, all nuclei reside in one of the two hyperfine ground 
states |<?i) and | g 2 ), with equal probability at room tem¬ 
perature. Further, we can assume that the nuclei of the 
different macroscopic ensembles l introduced above are 
evenly distributed among these two states as well. This 
collective ground state is denoted by |G). 

The definition of the collective excited states demands 
a more elaborate approach. In Ref. [19] collective ex¬ 
cited states E / j + ' 1 ) were introduced, which denote a sym¬ 
metrized excitation on the transition fi. Note that such 
an excitation is shared only by N/2 nuclei, since only half 
of the nuclei were originally in the ground state of the re¬ 
spective transition. Here, we now generalize these states 
and denote a collectively excited state in the ensemble l 
on the transition fi by More formally, we define 

it as 


H l} ) = 


1 


y/NW/2 


N {,} /2 

E - 


i4>(n) s™ I G) , (32) 


and again only half of the nuclei in the respective layer l 
contribute due to the ground state distribution. Each 
of the contributory nuclei couples to the modes with 
the same rate gMWl. This qualifies the collective states 
defined here to rewrite the system dynamics given in 
Eqs. (23)—(25) in the linear regime. In the process, the 
sums over the macroscopic number of atoms ^ is sim¬ 
plified to the much more manageable sum over the dif¬ 
ferent resonant layers l as 

E^” )[% S- } = E J\N^ c ^ ]{l} \EW){G\ . (33) 


This way, we obtain the effective equations for the linear 
regime 


H « = E ( 






2) (o y] c 


,blW 


x v /p^I|E«)(G|+H.c., 


R = Rc + Rn 


(34) 
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hls= y . £ (C-ix-i) (^‘V 111 ') 

At,*/ j,l,k 

x |ViVWJV{*> |-BW)(-B{ fc} | , (35) 

£cav[p] = X!E (C' 1 - 1 '^) 

A*,^ j,l,k 

X £[ P , |£?«)(G|, |G)(^ fc} |] . 

(36) 


Finally, the reflection coefficient reads 


Y^g m} *\/NW{E {l} \p\G) 


(®out‘®in) • 

(42) 


This set of equations will form the basis for the descrip¬ 
tion of the EIT experiment 1161 , which we will analyze in 

Sec. ivm 


C. Heuristic extensions 


R = Rc 


/ y / 2 

ain h v kU] + iA c ) 

x (o^ t -l x .d M ) ^/pVW (EW\p\G) . (37) 


B. Unmagnetized layers 


A commonly encountered scenario is the setting with¬ 
out a magnetic hyperfine splitting in the resonant layers. 
In this case, the level scheme of the 57 Fe nucleus reduces 
to a two-level system with one ground and one excited 
state. From the general theory above, this behavior can 
be emulated by setting the energy splittings 6 g , 6 e of 
the ground and excited states to zero and choosing the 
quantization axis such that the incident beam with polar¬ 
ization di n only drives the linearly polarized transitions 
At = 2,5 (c.f. Tab. 0- Further, we define the state 

I E {1} ) = -±= (| E\ l} ) + |£ 5 W >) , (38) 

which describes an excitation in the Zth resonant layer 
without the distinction of the two hyperfine substates. 
We obtain for the master equation and the reflection co¬ 
efficient 

Hq = + H.c. , 

j i 

(39) 

H h s = gU]{l}g[3 ]{ fe }* V / ivOh/VW 

3 l,k 

x \E {l} )(E {k} \ , (40) 

£ cav [p] = -E^iE-9 m{ V ]{fc} VivWivw 

j l,k 

x £[ P ,\eW)(G\,\G)(eW\] , (41) 


R = 


i+E 

j 


24 1 

Kb'] + i/± [ § 


i \f^ K ^R 

Oin “ K U] + iA [ £ 


Before we study the phenomena and consequences 
which emerge from our generalized theory, we take a step 
back and consider additional effects in our system, which 
in general will turn out to be of significance. We em¬ 
phasize, that so far our theory is developed to capture 
the cavity character of the layer system, i.e. the guided 
modes and the embedded resonant nuclei. However, in 
the grazing incidence geometry, also effects which are not 
related to the structure of the cavity and stem from bulk 
material properties become important. 

Since the refractive index of the cavity materials at x- 
ray energies is less than one, total reflection is observed 
for small incidence angles 9 in the few-mrad range, while 
for larger angles the light is essentially completely ab¬ 
sorbed. The transition between those two regimes is not 
sudden, but can be characterized by a smooth function 
7?Enveiope(#)- As soon as we consider the reflectance over 
a broader range of incidence angles, this envelope has to 
be taken into account. Note that previous studies have 
been performed at a fixed incidence angle mrnmm 
or covered only tiny angular ranges, for which the en¬ 
velope could be considered constant mi- To include 
this total reflection behavior at grazing incidence, we 
will heuristically combine the analytical formula from 
Eq. (301, describing the guided modes, with the reflection 
curve of the cavity’s mirror material only, which approx¬ 
imately takes into account the total reflection envelope. 
The envelope function i?Enveiope(^) is described in more 
detail in Appendix [Aj 

Furthermore, it has been observed that the cavity ma¬ 
terial dispersion leads to an additional relative phase be¬ 
tween light reflected off of the outside of the cavity and 
light entering the cavity. This dispersion phase was found 
to be necessary, e.g., to describe the asymmetry of the 
reflection curve R(9) around the minima of the guided 
modes [2B]. It can be included by generalizing the contri¬ 
bution which stems from the direct reflection on the cav¬ 
ity surface (“—1” in Eq. (30)) with an additional phase 
factor exp (iefic)- As a second heuristic extension, we in¬ 
clude such a phase phase factor as well, and allow for a 
complex variable r instead of the cavity surface ampli¬ 
tude — 1 with |r| ss 1. This has the additional advantage 
that such a modification can also take into account pos¬ 
sible effects of far off-resonant modes, which would give 
rise to a small constant offset to the reflection coefficient. 

With the heuristic modifications described above, the 




















cavity contribution to the reflection coefficients reads 


Rc(0) = -^Envelope (^) ( V + ^ 


2/C 


, (43) 


(44) 


/cM +iA^(0) 

with the cavity detuning (c.f. Eq. @» 

Ap (0) =w cos ( 9) 2 + sin (0^) - 1 

IV. EFFECT OF MULTIPLE MODES 


In this part, we will discuss the influence of multiple 
modes on the reflectivity, while we still restrict ourselves 
to a single thin layer of resonant nuclei in the cavity. 
Hence, the index l corresponding to the different nuclear 
ensembles in the cavity, will be omitted in the following. 

We start by considering the nuclear contribution to 
the reflectance only. Restricting ourselves to only one 
layer l, it can be easily seen from Eqs. (p4|)-(37l that 


the general form of contributions to the effective master 
equation does not depend on the number of cavity modes 
j. In particular, no new operators or additional couplings 
between the different collective states are present in the 
effective master equation, which is a direct consequence 
of the adiabatic elimination. The sole differences are the 
coefficients entering the expressions. For instance, gen¬ 
eralizing the driving Hamiltonian Hq from a single to 
multiple modes requires only the modification 


tog™ 


E 

3 






(45) 


on the level of a coefficient. Similar replacements of the 
variables are required for the other parts contributing to 
the master equation. 

With the knowledge that the basic equations in the 
cases with one and with multiple cavity modes j are 
equivalent, the results obtained in the original theory 
(c.f. Ref. [IS]) can be straightforwardly extended by re¬ 
placing the coefficients with their respective generalized 
counterparts. Doing so for the linear reflectance without 
magnetic hyperfine splitting and neglecting the trivial 
polarization dependency (a* ut -d; n ), this yields 


Riy = —i 


Since in the single-mode theory a Lorentzian profile 
was derived for the nuclear spectrum, we also recover 
this line shape here in the case of multiple cavity modes. 
Similar to the original model, it is shifted due to a collec¬ 
tive Lamb shift and broadened due to superradiance ma¬ 
lt can be seen from the denominator in Eq. (46) that each 



mode induces its own frequency shift and line broaden¬ 
ing. But typically for any angle of incidence, all but 
(at most) one mode are driven far off-resonantly as men¬ 
tioned above. Then, the according values for the cavity 
detuning Ac become large and their respective contribu¬ 
tions to the cooperative Lamb shift and to the superradi¬ 
ance diminish. From the numerator of the nuclear part, 
we find that the strength of the nuclear signal is typically 
determined by one dominant mode with the smallest Ac- 
Generally, the behavior is as follows. The nuclei can be 
excited by the external driving field via each mode j. 
This is represented by the first sum in the numerator 
of Eq. (461. The emission forms an independent second 


step and can again occur via each mode, as indicated by 
the second sum. Hence, in the general case, interferences 
between the different cavity modes can arise. Neverthe¬ 
less, the general Lorentzian structure of the line profile 
is unaffected by this and hence no qualitatively different 
features appear in the spectrum. 

The main difference to the single-mode result is found 
in the cavity contribution to the reflectance Rc given 
in Eq. (30), and accordingly its heuristic extension in 
Eq. @. Since we included multiple guided modes in 
the analysis above, it is clear that the resonances of these 
modes should become apparent in the reflection curve, 
i.e. when considering the cavity reflectance in dependence 
on the incidence angle 9. Indeed, the expressions we de¬ 
rived in Eqs. (30) and (43) highlight these resonances in 


its sum. The resonance of a guided mode j is encountered 
at 9 = 9q \ where = 0, and the reflection curve will 
exhibit a local minimum in the vicinity of the resonant 
angles. 

At this point it is of interest, how well the actual an¬ 
gular dependent reflection curve R(9) can be described 
by the cavity part of the reflectance. To this end, we nu¬ 
merically calculate the reflection curve using established 
semiclassical methods, such as Parratt’s formalism m 
or simulations by CONUSS mi, which both give equiva¬ 
lent results. As we will describe the EIT scenario from 
Ref- i.HI! below, we specialize to this particular cavity 
structure. The parameters of this cavity geometry are 
given in the middle column of Tab. [I] Note that the two 
resonant iron layers do not pose a challenge in this anal¬ 
ysis, since in this first step the nuclear resonances are 
omitted in the description of the angular dependent re¬ 
flection curve. 

In the specified cavity it is platinum that acts as 
cavity mirror material, therefore the envelope function 
-^"-Envelope(0) taking account for the total reflection is the 
reflectivity of a single infinitely thick Pt layer. We fitted 
the absolute value of Eq. (431 with a maximum number 


of cavity modes j = 5 to the expected reflection curve for 
the cavity, which was obtained using Parratt’s formalism. 
The fit was performed in the range 0 < 9 < 5 mrad and 
the obtained parameters are summarized in Appendix |B| 
The result is shown in Fig. [3j Clearly, the quantum opti¬ 
cal model together with the heuristic extensions are well 
suited to describe the reflection curve. We observe that 
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the last guided mode in Fig. [3] is not reproduced. How¬ 
ever, this is expected, since it is the sixth mode not in¬ 
cluded in the fit model, and since its angular position is 
not within the fit range. 

Interestingly, we find that also the phase behavior of 
the cavity reflection coefficient is reproduced very well. 
This is remarkable, because only absolute values were 
taken into account in the fit procedure. The analytic 
formula in Eq. (43) has only been corrected for a global 
phase to match the phase behavior predicted by the Par- 
ratt formalism. Considering the phases in more detail, 
a deviation can be seen at the second guided mode at 
9 sa 3 mrad. In contrast to Parratt’s formalism, the curve 
obtained with the quantum optical model features an ap¬ 
parent phase jump of 2ir at the resonance. To understand 
this artifact we note that in the vicinity around a guided 
mode, the reflection coefficient can be approximated as 
Rc ~ — 1 + 2k#/[k-MAc], which results in a minimum at 
the resonance angle where A c = 0. The reflectance van¬ 
ishes completely for 2k# = k, which is known as the crit¬ 
ical coupling condition. However, a residual reflectance 
occurs for both 2k# > k and 2k# < k, which correspond 
to the over- and undercritically coupled cases, respec¬ 
tively [55]. Looking solely at the modulus |I?c|> though, 
these two cases cannot be distinguished. The difference 
becomes only apparent when the phase of Rc around the 
resonance angle is considered: For an undercritically cou¬ 
pled cavity mode the phase remains in the same branch, 
however it undergoes an evolution to the next branch 
in the overcritically coupled case, which manifests as an 
apparent phase jump of 2-7r. Generally, it might be ben¬ 
eficial to not fit absolute values, but to use the complex 
values of the reflection curve instead. In this case also 
the over- and undercritically coupled modes should be 
captured correctly within the quantum optical descrip¬ 
tion. However, since we are interested mainly at the 
third guided mode later on, which is the mode at which 
the EIT spectra have been measured in Ref. [TB] , we will 
use the parameters obtained in the fit discussed above 
for our further analysis. 

Finally, it should be mentioned that it is not meaning¬ 
ful to extend the quantum optical descriptions to very 
large incident angles 9. On the one hand, the theoretical 
description of the perpendicular polarization directions 
might break down. On the other hand, distinct non¬ 
grazing incidence effects are expected, since the cavity is 
no longer probed in (000) Bragg geometry. The angu¬ 
lar cutoff should therefore be around the total reflection 
edge, which for typical cavity settings limits the number 
of guided modes to approximately five. 


V. EFFECT OF MULTIPLE RESONANT 
LAYERS 



FIG. 3. (Color online) Cavity reflectance as a function of 
the incidence angle 6. The minima denote resonance angles 
at which guided modes are driven resonantly. Including a 
dispersion phase and the heuristic extension of an envelope 
given by the topmost layer (gray dotted line), the quantum 
optical theory (blue, solid) can reproduce the exact result cal¬ 
culated with Parratt’s formalism (red, dashed) very well. In 
the quantum optical model only the first five guided modes 
were taken into account and the fit was restricted to the range 
0 < 9 < 5 mrad, as indicated by the shaded area. Although 
only the absolute value of the reflectance was fitted, the be¬ 
havior of the phases calculated with the two descriptions are 
very similar. 


index j in the coefficients throughout this section. 

Since the cavity reflection part Rc is unaffected by 
including multiple layers in the theory, we considering 
the nuclear contribution to the reflectance only. Starting 
from the general expressions given in Eqs. (34|-(37) and 
taking only one cavity mode j into account, we observe 
that the set of equations can be considerably simplified 
by means of a basis transformation. 

We introduce the states which resemble a collective 
excitation on transition /i, which is distributed among 
the different nuclear ensembles l, as 


In this part we will examine the influence of multiple ^ , - 

layers of resonant nuclei in the cavity, while restricting l^ + ^) = g E g^WN^lE^) (47) 

the analysis to only on cavity mode. Hence, we drop the i 
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with the normalization factor 

s= (Ek w r^ w ) 1/2 • (48) 

l 


Using the new states, Eqs. (34)—(37) become 


H n = (<C-U-Oin) |E< +} )(G| + H.c. , 

(49) 


Hls = ^6 hS g 2 J2(d^t ± -d^ c^c v \E^)(E^\ 


Eca.v \p\ = ~2^sS 'y \ i^dfi ‘ Iffi ' dij^j C^jCy 

fl.y 

x C[p,\eW)(G\,\G)(EW\] 
i t/kr 


(50) 


(51) 


R R ° a inK + i A c ^ C)1 


x (a: ut -l ± -d M ) (Ej+^\p\G). (52) 


Comparing these expressions with the effective master 
equation of the original single-layer theory (c.f. Ref. BSD. 
we observe an exact correspondence of the structure. 
Similar as in the case of the extension to multiple cavity 
modes, the differences manifest only in terms of the co¬ 
efficients. In particular, the collective coupling between 
the cavity mode and the nuclei is modified as 

\g\VN^g = (J 2 \ 9 {1} \ 2n{1} ) 1/2 > (53) 

i 


while all other relations remain the same. Hence, the 
shape of the reflection coefficient is unaffected by taking 
into account multiple resonant layers. In the absence of 
magnetization, it is given by 


R = 


-1 + 


2 kr 

k + iAc 


2 hr 


'(n + iAc) 2 A + q + lg^iCs-s^s) 


«. * ~ \ 

l out' a inJ • 

(54) 


VI. EFFECT OF MULTIPLE RESONANT 
LAYERS AND MULTIPLE MODES 


In the last Sections we generalized the theoretical de¬ 
scription to include multiple modes and multiple resonant 
layers, respectively. When restricting to one extension at 
a time, we observed that both give rise to additions in the 
nuclear reflection amplitude, while, however, leaving the 
general structure of a Lorentzian line shape unaffected, 
c.f. Eqs. (46 1 and (54). 


The general expressions for covering both multiple 
modes and multiple layers in the cavity were given in 
Eqs. (34) (37), or, in the absence of magnetization, in 
Eqs. (39) (42). To simplify this set of equations, it would 


be desirable to perform a basis transformation which con¬ 
verts the different states | E^) which describes an exci¬ 
tation in a single layer into a collective layer state, similar 
to the transformation we performed in Eq. (47). For that 


purpose, one would have to sum over the layers l, which 
then contains the coupling factor But since this 

coupling coefficient now also depends on the guided mode 
index j, the basis transformation must also involve the 
sum over the modes yHowever, it can be easily seen 


from Eqs. (39) (42) that this sum would be different for 


every contribution to the equations of motion, since the 
prefactors depending on j are mutually different. Hence, 
also in the absence of magnetization it is not possible to 
transform the system into a form in which only one col¬ 
lective state is excited. Rather, in a cavity configuration 
with l resonant layers the equations of motion need to 
be solved for the l coupled states | E^). This implies 
that the response of the nuclear ensemble will generally 
go beyond a Lorentzian line profile. 

The different coupling coefficients gbKT required for 
the extended theory need to be determined in different 
ways. Apart from a direct fit to numerical data, a more 
sophisticated approach is to derive the relative weights 
and phases from the field amplitudes calculated with Par- 
ratt’s formalism (2S1. This self-consistent method will be 


applied and explained in more detail in Sec. VIID 


VII. APPLICATION TO THE EIT SETTING 


Restricting ourselves to only one layer, the coefficient g 2 
reduces to | < 7 1 2 TV and we recover the result which we al¬ 
ready derived in Ref. US- 

Even though we included multiple layers in our analy¬ 
sis, we see from Eq. (54) that it is not possible to explain 
an EIT-like spectrum as reported in Ref. m- Rather, 
we will find that it is the combined extension of multiple 
layers and multiple guided modes to the theory, which 
will be able to explain the EIT phenomenon. This will 
be shown in the following Sections. 


Next, we will analyze a particular setting, in which 
both multiple layers and multiple modes are considered. 
In the previous Section, we already discovered that the 
general theory developed here allows for reflection spec¬ 
tra, which comprise nuclear responses beyond a simple 
Lorentz profile. As it was shown in Ref. m employing a 
semiclassical model, the reflectance of a cavity with two 
unmagnetized resonant layers can even exhibit EIT-like 
spectra. It is this particular scenario which will be dis¬ 
cussed in more detail in the following. 
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A. The EIT experiment 

The key of the setting studied in Ref. [16] is the place¬ 
ment of two ensembles of 57 Fe nuclei in the cavity. The 
incidence angle was chosen such that the third guided 
mode of the cavity is driven resonantly and the first layer 
of the resonant nuclei was placed in a field node, the sec¬ 
ond layer in a antinode as sketched in Fig. Ha). Fol¬ 
lowing the interpretation in the same reference, only the 
latter ensemble is probed by the x-ray beam and decays 
superradiantly with rate F, while the nuclei in the first 
layer are only subjected to natural decay 7 on a much 
longer timescale. However, the first ensemble can cru¬ 
cially influence the system’s dynamics, as a control field 
flc between the two layers is naturally established. It 
arises due to radiative coupling between the two ensem¬ 
bles. The resulting level scheme, as also visualized in 
Fig. 0a), is equivalent to a system featuring EIT. By 
employing a cavity like the one sketched in Fig. Ha), the 
key signature of EIT, transparency of the medium on 
resonance, could be verified in Ref. m- Interestingly, it 
was found that the coupling field f 1q vanishes by inter¬ 
changing the roles of the two layers, i.e. placing the first 
layer in a field antinode and the second ensemble in a 
field node, see Fig. H b )- 


B. Theoretical analysis 


In order to describe the setting of two resonant layers 
with our quantum optical model, we restrict ourselves to 
two ensembles of nuclei l = 1,2 in the cavity, but we 
still allow for an arbitrary number of cavity modes j. As 
before, we consider the linear response case without mag¬ 
netization and omit the trivial polarization dependence 
in the following. We rewrite the effective Hamiltonian 
from Eqs. (391 and (40) as well as the detuning part from 
Eq. (j 6 j) as 


H = + Ac.) + (h{ 2 >|e{ 2 >)(G| + Ac.) 

+ (<5 {1} - A)|E^ 1 } )(E {1} | + (<5 {2} - A)|A {2} )(£{ 2} | 


+ (£ { 1 ’ 2 } |£ { 1 } )(.E {2} | + Ac.) . (55) 


Here, the first line covers the driving of the two layers, the 
second line accounts for the cooperative Lamb shifts and 
the detuning, and the last line describes a coherent cou¬ 
pling between the two layers. Later, we will see that the 
last contribution can in parts be identified with the con¬ 
trol field Qc from the EIT interpretation in Ref. |T 6 ]. The 
incoherent Lindblad terms in our description are given by 

c = -(I + ^ <1} ) l^ {1} X G l’ l G ><^ {1} l] 

-(|+ 7 {2} ) £[p,\eW)(G\,\G)(eW\] 

- 7 {1 ’ 2} C[ P ,\eW)(G\,\G)(eW\] 

- 7 { 1 ’ 2 } *£[p, |E m )(G|, |G)(E {1} |] . (56) 


Here, the first line accounts for spontaneous emission and 
superradiance. The other two terms describe an incoher¬ 
ent cross-damping term H2130], which will contribute 
to the control field coupling in the EIT interpretation as 
well. The coefficients in Eqs. (55) and (56) are given by 


II 

-M 

<5 

( 57 ) 

m = Y / &g [jm \ 2 N {l} , 

( 58 ) 

5 {1 ’ 2} = £ ^ 2 gUmgUm^NWNm , 

( 59 ) 


( 60 ) 

yf 1,2 } = £ $ §^ 1{1 V 1{2}i Vivui Arm . 

3 

( 61 ) 


The effective level scheme of the system defined above 
is visualized in Fig. Hb). The similarity to the scheme 
used in the interpretation of Ref. m can already be an¬ 
ticipated. However, in our approach a larger number of 
coherent and incoherent couplings are present. Never¬ 
theless, the relative strength and hence the importance 
of the coupling rates can be straightforwardly estimated, 
as we will show in the following. 

As mentioned above, in the cavity geometries of inter¬ 
est, the 57 Fe layers are arranged such that one layer l = 1 
is located at a held node of the third guided mode, while 
a second layer l = 2 is located at an antinode. As a con¬ 
sequence, the nuclei in the node hardly couple to driven 
mode. In our quantum optical language, we can repre¬ 
sent this idealized case by setting the respective coupling 
constant to zero, i.e. gl 3 !! 1 ! = 0. At the same time, all 
other modes j / 3 are driven strongly off-resonant, such 
that their cavity detuning A^ becomes large. Indicating 
this suppression due to the large detuning with a sym¬ 
bolic notation 1/Ac, we find the scalings 


; AfU , ^{1,2} „ 

1 

Ac ’ 

(62) 


1 

A 2 ” ‘ 

^ c 

(63) 

In contrast, the coefficients 



^{2} , ^{2} , ^{21 

~ 1 

(64) 


for the second layer are not suppressed due to cavity 
mode detuning, as they still contain the non-zero cou¬ 
pling coefficient g[ 3 H 2 l to the resonantly driven mode. 
From these scalings we can already anticipate the EIT 
behavior in accordance with the interpretation discussed 
Ref. HU : Only the nuclei in the second layer decay super¬ 
radiantly. The collective decay of atoms in the first layer 
7 ^ and the cross-damping terms 7 O. 2 } are quadrati- 
cally suppressed in the detunings of the additional cavity 
modes and can be neglected in a first approximation. 
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FIG. 4. (Color online) (a) Effective level scheme of the EIT scenario from Ref. [16]. The nuclei in the field antinode are driven 
by the external field ftp and decay superradiantly with rate T, the nuclei in the field only due to spontaneous decay 7 . Both 
ensembles are coupled via the cavity field f ic, which gives rise to a scheme equivalent to EIT. (b) Effective level scheme in 
the quantum optical model. Two collective excited states are coupled to the ground state. Coherent driving and collective 
Lamb shifts are marked in blue, superradiant spontaneous emission is denoted by curly red single headed arrows, cross-damping 
between the excited states by curly red double headed arrows. [Thick solid / solid / dashed] lines denote the cavity mode 
detuning scalings ~ [1 / ^c 1 / ^c 2 ] and mark the relative magnitude of the different couplings. 


However, other contributions due to the presence of fur¬ 
ther cavity modes can have a substantial influence on 
the system, such as the coherent driving between 

the two layers, which can give rise to the coupling field 
required for EIT. 

These scalings with the cavity detuning Ac are visual¬ 
ized in the level scheme shown in Fig. [4j[b) as well. Cou¬ 
pling rates denoted by thick, solid or dashed lines and 
indicate the different powers in the scaling behavior with 
respect to Ac- With the relative magnitude of the rates 
in mind, a very close similarity with the EIT level scheme 
from Fig. |4|a) can be observed. Hence, our analysis so 
far also suggests EIT-like features in the system. How¬ 
ever, it is yet unclear how the additional driving terms 
and inter-layer coupling terms affect the spectrum in de¬ 
tail. In order to answer this question, we will now turn 
to the analytic solution of the model. 


Starting from Eqs. (55) and (56), we find that the equa¬ 


tions of motion for the density matrix elements 

Pig = (E {1} \p\G) , 

P2G = (E {2} \p\G) , 


(65) 

( 66 ) 


form a closed set of equations in the limit of linear re¬ 
sponse, i.e., where the populations (G|p|G) ~ 1 and 
(E^ IplE^ 1 }) = (E^\p\E^) S3 0 and the coherence 
between the excited states (E^\p\E^) vanishes. The 
equations of motion read 


dt 


dt 


Pig = 


i(A-5 {1} )-7 {1} -2 


Pig 


-inW-(i5W + jW)p 2G 


P2G ~ 


*( A-5 {2} )-7 {2} -3[ 


P2G 


-^ 2 >-(# 1 ’ 2 >*+7 {1 ’ 2} Vg- 


(67) 


( 68 ) 


From this we obtain the steady state solutions of the 
coherences 


Pig = 


P2G = 


A^m - (-ft 1 # + *7li>2}) 


AfHAf 2 ! -n 2 c 


AClHf 2 } - (-6IL2}* + ^{l,2}*) ft{l} 


AinAf 2 ! - n 2 c 


with the abbreviations 

A« =A-jW+iQ+7«), 


(69) 

(70) 

(71) 


n 2 c = 


(>’ 2 > - n {1 ’ 2} ) (^ {1 ’ 2} * - *7 {1 ’ 2} *) ■ (72) 


With the solutions for the coherences at hand, we can 
now turn to the observable, the complex reflection coef- 


2k 




ficient R. According to Eq. (421, it is given by 

* = -i + £- 

with 


K [j] + iA b] 


PlG + R {2} P2G , (73) 


# = --V 

(7:.. ^ 


2k 


[i] 


Atbl -p iA 


M 


2 g Um*y/N{i} . (74) 


At this point it is instructive to discuss the scaling related 
to the cavity detuning A <7 once again. As before, we as¬ 
sume that gl 3 !! 1 ! = 0, i.e. the first layer does not couple 
to the driven cavity mode j = 3 since it is located at a 
field node. In this case we find that ~ 1/Ac, while 
R^ 2 } is not suppressed due to a cavity detuning, since the 
second layer can couple to the resonantly driven mode as 
g[3]{2} ^ 0. Furthermore, from Eqs. (69) and (70) we find 
that pic ~ 1/Ac, whereas the pic is not suppressed in 
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this fashion. Therefore, for a qualitative understanding 
of the reflectance, it is well justified to drop the quadrati- 
cally suppressed contribution R\ 1 ^ pic and only consider 
the reflection signal which stems from the second layer, 
i.e. R^p 2 G- 

In a further step, we restrict the numerator of the frac¬ 
tion in P 2 G to terms up to linear order in 1/A c- 
Moreover, we neglect the tiny collective Lamb shift and 
superradiance of the nuclei in the first layer. This yields 
the reflection coefficient 


* = -i + £ 

3 


24 ' 

rebl + 


+ r { 2} n 2 


_ A + ij 

(A + i 2) (a - S 2 + * + 72)) 


(75) 


The nuclear contribution to the reflectance is revealed in 
the second line. Its spectral shape is essentially that of a 
system featuring EIT. Hence, we recover the same result 
as in Ref. [16] : In a cavity with two resonant layers it is 
possible to realize the phenomenon of electromagnetically 
induced transparency. 


C. Comparison to the semiclassical analysis 


In Ref. 116) . the case was studied, where the empty- 
cavity contribution to the reflectance vanishes, and a 
semiclassical theory based on transfer matrix techniques 
was used to derive an expression for the nuclear re¬ 
flectance. In the sign convention of the nuclear reso¬ 
nances of the present work, the result was found as 


R — -*^2/02^2—b 


a + q 


(A + *i)(A + i£)-ft£ ’ 


r = 7(1 + d2/o£?2—) , 

Y = did2fo\E 2 -+E 1+ _ , 


(76) 

(77) 

(78) 


where d\ and d 2 are the thicknesses of the respective 
two layers, /o is the nuclear scattering amplitude at reso¬ 
nance, and E 2 _f, E 2 _, and _are transfer matrix el¬ 

ements. Comparing it with the part of the nuclear reflec¬ 
tion in the quantum optical expression given in Eq. (751, 
we notice a perfect agreement of the structures of the two 
formulas. However, as an important consistency check, 
it remains to be verified if the scaling with the number 
of nuclei in the two layers agrees as well. In the semi¬ 
classical theory it was shown that the amplitude of the 
reflection coefficient and the superradiance of the nuclei 
in the second layer scale linearly with the thickness of the 
second layer d 2 , and furthermore the control field Qc was 
shown to be proportional to y/did 2 - The present model 
does not directly contain the layer thicknesses as param¬ 
eters. But since d\ oc N W and d 2 oc N^- 2 \ is sufficient to 
show that the scaling relations also hold for the numbers 


of nuclei. From Eqs. (57) (61) and (74) it can indeed be 


seen that the relations are correctly reproduced by our 
theory. 

This is an important result, since it is a strong hint 
that the two independently derived results do not coin¬ 
cide by chance, but also agree on a more fundamental 
level. Hence, the model developed here can now be em¬ 
ployed to shine light on the EIT scenario from a different 
perspective. 


In the nuclear reflectance calculated in Eq. (75), the 
coupling Rabi frequency occurs as in the denomina¬ 
tor, whereas in standard EIT settings it appears as a pos¬ 
itive real-valued variable |flc| 2 - Taking a closer look at 
our definition of the coupling Rabi frequency in Eq. (72), 


we note that in our case 12^ can generally be complex. 
Also in the semiclassical theory the complex field ampli¬ 
tudes and transfer matrix elements E 2 _ \.E\^ _allow for 

complex values, c.f. Eq. (78). The results of Ref. m, 
though, seem to imply that the imaginary component is 
very small and an EIT situation is well realized. How¬ 
ever, from the theoretical analysis of the semiclassical 
models, this fact could not be understood and the influ¬ 
ence of the imaginary component was unclear |25| . With 
the present theory, though, it is now possible to exam¬ 
ine the complex nature of the coupling in more detail. 
From Eq. (72) we know that it is not only given by the 


coherent coupling between the two layers as writ¬ 

ten in the Hamiltonian in Eq. (55), but is also affected 
by the incoherent cross-damping term yl 1 ’ 2 ! between the 
two layers. In the discussion on the scalings we have al¬ 
ready seen that, in contrast to the coherent contribution, 
the incoherent term is suppressed quadratically with the 
detuning of the off-resonant cavity modes. Thus, we find 
that the incoherent part 


Im(II^) = —2 Re 


}* 


W (79) 


can be neglected, and ~ |^ 1,2 ^| 2 such that the real 
component of the coupling frequency flc dominates. 

Furthermore, the microscopic ansatz of our quantum 
optical theory enables one to interpret the origin of the 
coupling between the layers. While in Ref. m it was 
shown that the EIT control field arises from radiative 
coupling between the two resonant layers, it can now be 
pinned down from Eqs. (59), (61) and (72) to 


n 2 c = (%) 2 nWnW 

' gUmgbm* 


y a 



g vm*g\jm 
Ag ] - mb] 


(80) 


Since we assumed that the first layer does not couple to 
the third guided mode in the idealized case, i.e. t/ 3 !! 1 ! = 
0, we observe that the coupling field is only mediated via 
the remaining guided modes j / 3 in the cavity. This 
way, it becomes now also clear why the EIT phenomenon 




















14 


was not obtained in Sec. [Vj where multiple layers, but 
only one guided mode was included in the theoretical 
analysis. 


D. Numerical analysis 


Let us now see how well our analytical expression for 
the reflectance derived above performs in practice. In 
particular, we aim to describe the spectrum of the EIT 
cavity defined in the first column of Tab. [T] with our quan¬ 
tum optical model. Moreover, we include a second cavity 
into the analysis: While the EIT cavity has its resonant 
layers in a node and antinode of the field of the resonantly 
driven mode, respectively, we also consider a cavity in 
which the situation is reversed. Namely, the first reso¬ 
nant layer is located at a field antinode and the second 
ensemble of nuclei at the field node. The corresponding 
geometry is defined in the last column of Tab. [T] The 
two cavity layouts reflect the cases discussed in Ref. M, 
where it was shown that the first cavity exhibits the EIT 
phenomenon, while for the second system the control cou¬ 
pling £lc vanishes and only a Lorentz-like spectrum is 
measured. 

In order to determine the free parameters related to 
the cavities defined in Tab. |T] for the quantum optical 
model in a consistent way, we employed the following 
method. First, we restricted ourselves to the first five 
guided modes in the theory and did not take into account 
the resonant nuclei yet. For each of these modes the 
angles 0$, at which the modes are driven resonantly, 
and the decay and coupling rates and have to 
be determined. The parameters can be found by fitting 
Eq. (43) to the reflection curve as function of the x-ray 


incidence angle as it was already done in Sec. IV 


With the cavity parameters at hand, the next step 
now is to include the nuclear resonances to the model. 
In particular, the complex collective coupling coefficients 
gtiW} v / 7Vf z 1 between the jth guided mode and the layer l 
of resonant nuclei have to be determined. For each cavity 
the number of coupling coefficients is 10, since we special¬ 
ize to five guided modes in the analysis and each mode 
can couple to the two respective layers. In order to avoid 
arbitrariness in a fit to numerical data, it is advisable to 
decrease this large number of free parameters. In fact, is 
it possible to determine all coupling coefficients in a con¬ 
sistent way, while keeping only one global scaling as free 
parameter. To illustrate this, we note that the coupling 
coefficients can be decomposed as 


gUm^/NW = £\JlW • (<7«A/iV«) , (81) 

where the first factor denotes the cavity field amplitude 
of mode j at layer l, and the second factor includes the 
collective nuclear dipole moment. Next, we exploit that 
the complex field amplitudes in the cavity can be eas¬ 
ily derived by means of Parratt’s formalism [25]- In a 
simple picture, we can interpret the resonant nuclei in 


the cavity as a perturbation, which modifies the cavity 
field and, accordingly, the reflectance. The cavity field 
in the presence of nuclear resonances can be understood 
as a superposition of the bare cavity field and the con¬ 
tribution due to scattering at the nuclei. This presuppo¬ 
sition clearly holds for x-ray frequencies apart from the 
nuclear resonance, but does also give consistent results 
directly at the resonance, where the perturbation due to 
the nuclei is not generally small. Hence, to determine the 
field coefficients , it is not necessary to include any 

nuclear resonances in Parratt’s formalism, but only the 
bare cavity field in the absence of 57 Fe resonances are 
required. With the input field normalized to intensity 
one, like in Fig. [2] we can directly compute all complex 
valued field coefficients £b']P} at the center of the respec¬ 
tive layers l by tuning the incidence angle 9 to the angles 
9^\ where the jth cavity mode is driven resonantly. The 
remaining task is to determine the second coefficient in 
Eq. ( [8l] ), which takes into account the nuclear proper¬ 
ties and other constant contributions. Since both iron 
layers in the cavities have the same thickness and hence 
the number of nuclei is the same, we can expect that 
gWy/NW is a constant and acts only as a scaling pa¬ 
rameter for the previously determined field amplitudes. 
This way, all coupling coefficients (/b']W-\/Ar{T can be de¬ 
duced by fitting the model to numerical data, calculated 
with Parratt’s formalism, with only one free scaling pa¬ 
rameter. The complex field amplitudes for both cavities 
and the scaling parameter found in our analysis are sum¬ 
marized in Appendix [B] We note that the couplings to 
the layers, which are located in the field nodes of the third 
cavity mode, do not completely vanish due to the finite 
thickness of the layers and a potential misplacement in 
the cavity. However, they are found to be much smaller 
than the coupling coefficients of the respective layers in 
the cavity field antinode. This can already be deduced 
from the field intensity distributions shown in Fig. [2j 
An alternative approach to determine the coupling 
constants gblW\/iVW is by fitting the model with all 
coefficients directly to numerical data. While this proce¬ 
dure is not as persuasive as the consistent method de¬ 
scribed above, it might also offer some advantages in 
quantitative studies. Errors in other parameters, such 
as the coefficient which characterizes the cavity modes 
can partly be compensated. Moreover, for iron layers 
with a larger thickness the field amplitude might not be 
constant and, in contrast to the method from above, an 
effective coupling strength would be naturally obtained. 
Finally, the fitted parameters could provide a handle to 
cover the fact in more detail, that on-resonance the nuclei 
have an effect on the cavity field which goes beyond a per¬ 
turbation. In this work, however, we will not optimize 
the parameters in this way, but utilize the coefficients 
derived previously to illustrate the general consistency of 
our model. 


Now we are able to benchmark our analytical result for 
the case of two resonant layers, which was calculated in 
Eqs. (69) (74). A comparison with the frequency- and 
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FIG. 5. (Color online) The reflectance of (a) the EIT and (b) the non-EIT scenario is shown as a function of the detuning A and 
the incidence angle 6. The results derived with the extended quantum optical model agree very well with the predictions from 
Parratt’s formalism. The dashed line at 6 as 3.5 mrad marks the angle at which the 3rd cavity minimum is expected. A cut 
along this line corresponds to the spectrum measured in Ref. [16] and is shown in Fig. [6] Parameters are given in Appendix |B| 


angular-dependent reflectance for the EIT and the non- 
EIT cavity is shown in Fig. [5] Clearly, the agreement 
between the two different models is very good. We stress 
that this is not an obvious result, since the parameters 
for the quantum optical model were determined indepen¬ 
dently and not obtained from a fit to the numerical data. 

A range, in which strong deviations can be observed, 
is the domain around A « 0. Here, the exact numeri¬ 
cal solution obtained from Parratt’s formalism shows an 
additional structure. This can be understood from the 
following considerations. If the x rays are not resonant 
to the transition in the 57 Fe nuclei, they will primarily 
be damped due to the electronic absorption in the cav¬ 
ity, before they can reach the lower resonant layer. If, 
however, their frequency is too close to resonance, the x 
rays will additionally be absorbed by the nuclei in the 
upper layer. Consequently, the field seen by the nuclei 
in the second layer is strongly modified compared to the 
off-resonant case. However, in the derivation above we 
assumed that the presence of the nuclei can be treated as 
a small perturbation to the cavity field, which is not the 
case in the extreme situation encountered here. An ap¬ 
proach for future studies could thus be to comprise this 
effect self-consistently into the quantum optical theory 


for an even better agreement with the numerical data. 

We now turn to the spectrum measured at the inci¬ 
dence angle corresponding to the third guided mode, i.e., 
the situation from Ref. m- The spectra for both the 
EIT and the non-EIT cavity defined in Tab. [I] are shown 
in Fig. [6] Again, we observe a good qualitative agree¬ 
ment of our theory with the numerical data obtained with 
Parratt’s formalism, which could already be anticipated 
from the accordance in Fig. [5] But in any case, the fact 
that the EIT as well as the non-EIT spectrum is repro¬ 
duced without post-optimization of the consistently de¬ 
rived parameters, supports the validity of our theoretical 
description. 

Finally, we want to review to role of the coupling field 
Hp. In the theoretical analysis in Sec. |VIICl it was found 
that the presence of this control field gives rise to EIT. 
Moreover, we found that the control field is established 
by an interaction between the nuclear ensembles via dif¬ 
ferent cavity modes, as the layer in the field node does 
not directly couple to the driven guided mode. From 
our numerical analysis we observe that this idealized case 
is not strictly realized. Since the coupling coefficient is 
small yet finite, also the resonantly driven mode gives 
rise a coupling between the two layers in the cavity. Fur- 
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FIG. 6 . (Color online) Spectra of (a) the EIT and (b) the 
non-EIT cavity at incidence angle 9 ss 3.5 mrad, at which 
the third guided mode is excited. The quantum optical de¬ 
scription (blue solid lines) is in qualitative agreement with the 
exact result derived with the Parratt formalism (red dashed 
lines). Parameters are given in Appendix [5] 

thermore, the control field Vic does not vanish in the 
non-EIT case and hence the Lorentz-like spectrum can¬ 
not be explained by its absence in the frame of our model. 
Rather, in the non-idealized case it is the interplay with 
other contributions to the reflection coefficient and their 
interference which results in the Lorentzian spectrum. 


SUMMARY AND DISCUSSION 

In summary, we investigated the effect of multiple 
modes and multiple ensembles of resonant Mossbauer nu¬ 
clei in an x-ray cavity QED setup, which has recently 
served as a platform for multiple experiments related to 
x-ray quantum optics. Most of the time, the scenario 
with a singe ensemble, realized by a layer of collectively 
acting nuclei, has been studied theoretically jj9] as well 
as in several experiments [10HT3] . The theoretical frame¬ 
work of this work is applicable to also model experimental 
settings with more than one resonant ensemble m and 


interpret them from a quantum optical point of view. 

Our theory from Sec.[TTT]is based on the approach taken 
in Ref. }l9j and constitutes a generalization with multiple 
cavity modes and several layers of resonant Mossbauer 
nuclei. Similar to the original theory, we were able to 
simplify the basic equations using two well justified ap¬ 
proximations. By adiabatically eliminating the cavity 
modes and restricting the analysis to the linear regime, 
effective equations of motion for the nuclear ensembles 
could be derived. The resulting set of equations charac¬ 
terizing the dynamics of a few-level system can easily be 
solved analytically. 

In Sections [TV] and |V| we discussed the consequences of 
the two extensions to the theory in detail. By introduc¬ 
ing multiple cavity modes to the model we found that the 
spectral properties around the resonance of 57 Fe are un¬ 
affected. In absence of magnetic hyperfine splitting, the 
nuclear response is given by Lorentz profiles, which are 
shifted and broadened due to collective effects. The dif¬ 
ferences to the predictions from Ref. [T5] manifest only 
in the coefficients entering the final expressions. How¬ 
ever, a clear difference could be observed when the re¬ 
flectance was studied as a function of the x-ray incidence 
angle. While a single-mode theory can only indicate one 
guided mode of the system at a time, our extension al¬ 
lows to accurately model the reflectance over a range of 
several mrad, reproducing all guided modes. Moreover, 
we found that our model, which takes into account the 
effect of the cavity and its modes, can be heuristically 
extended to incorporate bulk material properties such as 
the total reflection envelope. This way, a close agreement 
with established semiclassical models could be achieved. 
Further, the effect of multiple ensembles of 57 Fe nuclei in 
the cavity, located in different layers, was studied. This 
extension alone did not give rise to qualitatively new ef¬ 
fects. 

Next, we analyzed the case in which both extensions 
enter the theory at the same time, i.e. multiple cavity 
modes and multiple resonant layers. We could show that 
in this case the equations cannot be mapped to an effec¬ 
tive two-level system, as the coupling coefficients between 
the different nuclear ensembles and cavity modes are mu¬ 
tually different and do not allow for a diagonalization in 
which only one excited state is probed. Rather, more 
advanced level schemes generally occur in this setting. 

In the final part of this work we applied the general 
theory to the setting, which was experimentally explored 
in Ref. [16]. In this reference, EIT-like spectra could be 
observed for a cavity with two layers of resonant iron 
nuclei. We applied our quantum theoretical approach to 
the setting and could successfully reproduce the findings. 
An effective level scheme with one collective ground and 
two collective excited states which captures the complete 
system could be found and an analytic solution for the 
reflection coefficient was given. For the idealized case 
of perfectly placed layers in the cavity we found that 
the nuclear response has indeed the spectral shape of a 
system featuring EIT. In this process, we compared our 








17 


result to the previously used semiclassical models and 
observed agreement on the analytic level. In particular, 
the scalings with respect to the number of atoms in the 
respective ensembles are reproduced by our quantum op¬ 
tical description. Most importantly, the question on the 
nature of the control field, which forms a pivotal require¬ 
ment of EIT, could be elucidated. From our analysis of 
the idealized scenario we found that the radiative cou¬ 
pling is mediated by the off-resonant cavity modes. 

We further developed a way to consistently derive the 
different coupling rates required for the model. This ap¬ 
proach is based on an analysis of the cavity in the ab¬ 
sence of nuclear resonances. Hence, the arising spectral 
features can be traced back to the capability of our model 
and are not due to a potentially biased parameter fit. In 
our numerical data we observe a good agreement to the 
results of semiclassical models and the essential features, 
such as the signatures of EIT, are reproduced. 

While we mainly analyzed the cavity properties in the 
absence of magnetization, we emphasize that the ex¬ 
tended theory description developed in this work is not 
restricted to a vanishing magnetic hyperfine splitting in 
the resonant layers. Rather, in our model it is possi¬ 
ble to include all Zeeman sublevels properly. In future 
works, this could be exploited to combine the effect of 
multiple layers and modes, giving rise to the EIT-like 
effects, and magnetization, leading to the phenomenon 
of spontaneously generated coherences PE]. This way, a 
broad class of quantum optical level schemes could be en¬ 
gineered, indicating promising perspectives of x-ray cav¬ 
ity QED with Mossbauer nuclei. 


man National Academic Foundation. 


Appendix A: Reflection curve envelope due to bulk 
material properties 

As explained in the main text, the quantum optical 
theory can model the minima in the reflection curve R(0) 
indicating the guided modes of the cavity, whereas the 
envelope formed by the total reflection behavior is not 
part of the description (gray dotted line in Fig. [3]). How¬ 
ever, it can be included by combining the expressions of 
the quantum optical model, describing the cavity struc¬ 
ture, with an envelope function i?Enveiope(0), which takes 
into account the reflection of the bulk material. This 
way, a good agreement to the semiclassical approaches is 
achieved. For a single (infinitely thick) layer, the Fresnel 
reflection coefficient reads 

sin (0) — \/sin (0) 2 + n 2 — 1 

^Envelope (6) = - , > (Al) 

sin (0) + yj sin (0) 2 + n 2 — 1 

where 0 is the angle of incidence and n is the refractive 
index of the material. For a platinum layer and at x-ray 
energy 14.4 keV, the latter is given by p|T 

n = 1 — 8 + i/3 , (A2) 

S = 1.603365 x 1(T 5 , (A3) 

f3 = 2.56353 x 1CT 6 . (A4) 
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In Tabs. IV and [V] we summarize the parameters for 
our quantum optical theory, which we used in this work. 


TABLE IV. Parameters for the EIT cavity. 


Mode j 

9q ' [mrad] 

[ 7 ] 

4 ! 

£\j]{ 1 } 

£[j}{ 2} 

1 

2.55943 

145807 

6110 

0.609 + 0.036 i 

0.903 + 0.346 i 

2 

2.99211 

533322 

311376 

2.105 + 1.795 i 

0.818 + 0.810 i 

3 

3.54936 

615909 

275736 

—0.031 + 0.440 i 

-1.683- 1.815 i 

4 

4.14850 

783648 

373532 

-0.947- 1.130 i 

—0.244 — 0.781 i 

5 

5.07939 

1718031 

767833 

-0.370 - 1.567 i 

0.361 + 1.582 i 


Note that the coefficients £ PH 1 ! for the EIT cavity and 
£[3]{ 2 } f or |q ie 

non-EIT cavity are near zero, indicative 
of a field node at the location of the respective layers. 
The asymmetry parameters r for the two respective cav¬ 


ities, expected to be close to the value — 1, were deter¬ 
mined as 7~eit ~ —0.981 + 0.363 i ss —1.047 e -0 - 354 * and 
^non-EIT ~ -0.979 + 0.383 i « -1.051 e ~ 0 - 373i . For the 
scaling parameter of the coupling coefficients the values 
pUl/vU} = gWy/NV} = 1983.897 were used. 
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TABLE V. Parameters for the non-EIT cavity. 


Mode j 

6q^ [mrad] 

[ 7 ] 

4 ] [7] 

gli]{ 1 } 


1 

2.58446 

242554 

150696 

1.889+ 1.613 i 

0.957 + 0.935 * 

2 

2.96021 

311202 

48680 

-0.760 - 0.140 i 

-1.648- 1.160 i 

3 

3.55108 

607732 

342667 

-1.631 - 2.086 i 

—0.058 — 0.332 i 

4 

4.17107 

922367 

453227 

0.203 - 0.067i 

1.031 + 1.606* 

5 

5.09251 

1844998 

833796 

0.217+ 1.585* 

-0.203- 1.604* 
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